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This is the third of a series of papers where we in- 
vestigate the noisy Burgers equation in the context of 
modelling a growing interface; for a brief account of the 
present work we also we refer to [Q . In the previous two 
papers, in the following denoted paper I U and paper II 
j|; a brief account of paper II also appeared in M, we 
discussed the originally proposed one dimensional noise- 
less Burgers equation ^ (J from a solitonic point of view 
and the noisy one dimensional Burgers equation [f7|8| in 
terms of a Martin-Siggia-Rose path integral pf-p], re- 
spectively. 

Phenomena far from equilibrium are widespread, in- 
cluding turbulence in fluids, interface and growth prob- 
lems, chemical reactions, biological systems, and even 
economical and sociological structures. In recent years 
much of the focus of modern statistical physics and 
soft condensed matter has shifted towards such systems. 
Drawing on the case of static and dynamic critical phe- 
nomena in and close to equilibrium |l5|]l6| , where scaling, 
critical exponents, and universality have served to orga- 
nize our understanding and to provide calculational tools, 
a similar approach has been advanced towards nonequi- 
librium phenomena with the purpose of elucidating scal- 
ing properties and more generally the morphology or pat- 
tern formation in a driven state. 

In this context the noisy Burgers equation provides 
maybe the simplest continuum description of an open 
driven nonlinear system exhibiting scaling and pattern 
formation. The Burgers equation was originally sug- 
gested in the one dimensional noiseless version [^|,[||l7] [p]] 



as a model for irrotational hydrodynamical fluid flow. 
Here u is the irrotational velocity field and v a damping 
constant or viscosity. Choosing the nonlinear coupling 
A = — 1, we recognize the usual nonlinear convective 
term p0| . As a nonlinear model for hydrodynamical tur- 
bulence [ET1- 
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has been studied intensively. 
It was early recognized that the nonlinear structure of 
the damped velocity field is dominated by shock waves 
|^|,^9|, yielding an inverse cascade and that the Burg- 
ers equation does in fact not characterize Navier-Stokes 
turbulence which is governed by a direct cascade and Kol- 
mogoroff scaling in the inertial regime. 

Treating the field u as the local slope of a growing 
interface, we analyzed in paper I, where a more com- 
plete bibliography can be found, the Burgers equation 
from the point of view of a soliton-carrying nonlinear 
damped evolution equation [ p8| . The equation describes 
the transient damped evolution of an initial slope configu- 
ration uq in terms of a gas of moving right hand viscosity- 
smoothed solitons connected by ramp solutions and with 
superposed damped diffusive modes with a gap in their 
spectrum proportional to the soliton amplitudes. In a 
heuristic sense the transient morphology is thus char- 
acterized by two kinds of excitations: Nonlinear soliton 
modes and linear gapful diffusive modes. This picture is 
also borne out by the nonlinear Cole-Hopf transformation 
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relating the slope field u to the diffusive field w satisfying 
the linear diffusion equation 



dw 
~dt 



with solution 



W (x,t)=Jdx' G(x - x',t)w (x') 
G(x,t) = [47wtp5 exp [-x 2 /Avt] , 



wq(x)= exp 
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allowing for a steepest descent analysis in the inviscid 
limit v — ► HH§. 

The emphasis on the nonlinear modes in paper I and 
the insight gained formed the natural starting point for 
our study of the noisy Burgers equation in paper II. This 
equation has the form ||7|,p|pl|-|33]] 



^- - AuV I u = vV z u + V?i , 
at 
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where the Gaussian white noise driving the equation is 
spatially short range correlated according to 



(rj(x, t)r)(x',t')) = A6{x - x')6(t - t') , 



(1.8) 



characterized by the noise strength A. Equation (1.7) has 
a much richer structure than the noiseless counterpart. 
The noise is here a singular perturbation in the sense that 
even a weak noise strength A e ventually drives the mor- 
phology described by Eq. (1.7) into a stationary driven 
state; the characteristic time scale is set by t oc ln(l/A) 
which diverges for A — > 0. The singular nature of the 
noise is also reflected in the known stationary probabil- 
ity distribution for the slope field |3l|] 



P st (u) oc exp 
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Here A constitutes an essential singularity in the same 
manner as the temperature entering the usual Boltzmann 
factor. 



Recasting the stochastic Langevin equation (1.7) in 
terms of a Martin-Siggia-Rose path integral [p|- |12| , [14| , [l3| 
we proposed a principle of least action in the nonper- 
turbative weak noise limit A — ► and derived canonical 
saddle point or field equations (note that in paper II the 
noise field ip was rotated, ip — > — itp) 
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coupling the slope field u to a deterministic aux iliary 
noise field ip. The coupled field equations (1.1C) and 



(1.11) effectively replace the stochastic equation (1.7) and 
describe the morphology of a growing interface. In ad- 
dition to the right hand soliton already present in the 
noisel ess ca se as discussed in paper I, the field equations 
(1.10) and (1.11) also admit an equivalent left hand soli- 
ton. A growing interface can thus be viewed as a gas 
of connected right hand and left hand solitons with su- 
perposed diffusive modes. The weak noise approach also 
allows a dynamical description of the soliton and diffusive 
mode configurations and associates an action S, energy 
E, and momentum II with a particular interface morphol- 
ogy. The statistical weight of a configuration is given by 
exp [— Sy A] with the dynamical action S, providing the 
generalization of the Boltzmann factor exp [— E/T] for 
equilibrium processes to dynamical processes. 

The Martin-Siggia-Rose path integral approach more- 
over permitted a simple interpretation of the scaling 
properties of a growing interface. The dynamical expo- 
nent z — 3/2 , entering in the dynamical scaling relation 



(u(x,t)u(0,0)) = \x\^- 2 F(\t\/\xn , (1.12) 

is thus related to the gapless soliton dispersion law 

E oc AIT , (1.13) 

whereas the roughness exponent £ = 1/2 follows from the 
spectral representation 



(u(x,t)u(0,0)) 



dK G(K) exp [-iEt + iKx] 
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assuming that the form factor G(K) ~ const, in the scal- 
ing region K ~ 0. The dynamical scaling universality 
class is associated with the lowest gapless excitation, i.e., 
for A ^ the soliton mode. For A = the Edwards- 
Wilkinson (EW) universality class |3l|,|3(],|lJ emerges 
with a gapless diffusive mode u> — vW , corresponding 
to z = 2, Q = 1/2 being unaltered. Furthermore, we 
derived a heuristic expression for the scaling function F 
in terms of the probability distribution for Levy flights 
ifHIJIi?! ; the scaling function has also been accessed by a 
mode coupling approach 

Summarizing, the weak noise saddle point approach to 
the noisy Burgers equation advanced in paper II yields a 
many body description of the morphology of a growing in- 
terface in terms of two kinds of quasi-particles or elemen- 
tary excitations: Nonlinear soliton modes corresponding 
to the faceted steplike growth of an interface with super- 
posed linear diffusive modes. Furthermore, the scaling 
properties and the notion of universality classes follow as 
a byproduct from the dispersion law of the lowest gapless 
excitation. For details and references we refer the reader 
to the somewhat tutorial presentation in paper II. 

Whereas a good understanding of the one dimensional 
case has been achieved both by renormalization group 
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methods Qj§,^l|-P3Lg7| and |4j^,|4J, by mappingto di- 
rected polymers |31[| 7by mapping to spin chains [|48|-p0| , 
from the lattice exclusion model f[5l|-[53T] , and by the soli- 
ton approach in paper II, the general case in d dimensions 
has proven much more difficult. In d dimensions the noisy 
Burgers equation takes the form fii\M 
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Here the longitudinal vector field u n , n = l,..d, is as- 
sociated with the height profile h of a growing interface 
according to (l n is a line element) 



U n = V„ft , 
h(x n )= / dl n U n ■ 



(1.16) 



Furthermore, V 2 = V„V„ and we assume summation 
over repeated Cartesian indices. The height field h is 
thus the underlying pote ntial for th e for ce or slope field 
u n - It follows from Eqs. ( |l.l5| ) and ( 1.16 ) that h satisfies 
the Kardar-Parisi-Zhang (KPZ) equation ]3l|-|33|] 



— = vV 2 h + -V n hV n h 
dt 2 



F 



(1.17) 



for a growing interface in d dimensions. Assuming (r/) = 
we have for completion here introduced the drift term 
F = -(\/2){V n KV n h) in Eq. ( |1.17| ) in order to en- 
sure that (h) decays in time in a co-moving frame. In 
Eqs. ( |TTT5| ) and ( ]L17| ) v is damping constant or viscosity 
characterizing the linear diffusive term, A is a coupling 
strength for the nonlinear mode coupling or growth term, 
and, finally, r\ is a Gaussian white noise driving the equa- 
tion and correlated according to 

(r,(x n ,t)r)(x' n ,t')) = Al[6(x n -x' n )5(t-t') , (1.18) 



where A is the noise strength. 

In higher dimension dynamic renormalization group 
calculations ||l] |33|] yield a (lower) critical dimension at 
d = 2 and a kinetic phase transition above d = 2, sep- 
arating a smooth phase characterized by the EW uni- 
versality class yielding Q = (2 — d)/2 and z = 2 and 
a rough phase characterized by a strong coupling fixed 
point. On the transition line renormalization group cal- 
culations and scaling arguments based on the mapping to 
directed polymers yield the exponents £ = and z = 2 
and suggest an upper critical dimension d = 4 jHJ . Most 
recently, an operator expansion method has been applied 
to the strong coupling phase yielding ((, z ) — (2 /5, 8/5) 
in d = 2 and ((, z) = (2/7, 12/7) in d = 3 g|§. 

In the present paper we focus on the stationary and 
time-dependent probability distributions for the height 
and slope fields described by the Burgers and KPZ equa- 
tions (1.15) and (1.17), respectively. As discussed in 
paper II these distributions are basically given by the 



Martin-Siggia-Rose path integral weighted by the effec- 
tive action for the appropriate paths. In the weak noise 
limit only the paths governed by the saddle point equa- 
tions contribute to the distributions and, as will be dis- 
cussed here, we can actually circumvent the path integral 
formulation entirely by a more direct approach based on 
the Fokker-Planck equation. In the weak noise limit this 
equation takes the form of a Hamilton-Jacobi equation 
implying a symplectic structure and immediately lend- 
ing itself to a canonical phase space formulation. We are 
thus able in a very direct manner to map the stochas- 
tic processes described by the KPZ-Burgers equations 
to a conserved dynamical system with orbits satisfying 
deterministic canonical Hamilton equations, identical to 
the saddle point equations in the path integral approach. 
The stochastic nature of the Langevin equations is re- 
flected in the peculiar topology of the energy surfaces. It 
turns out that the stationary probability distribution is 
determined by an infinite-time orbit on the zero-energy 
manifold whereas the time-dependent distribution, ap- 
proaching the stationary one at long times, corresponds 
to a finite-time orbit. Below we highlight some of our 
results. 

(i) In the generic case of a general nonlinear Langevin 
equation for a set of stochastic variables driven by white 
noise, the weak noise limit of the associated Fokker- 
Planck equation takes the form of a Hamilton-Jacobi 
equation, which in turn implies a symplectic structure 
with a principle of least action, an action, an associated 
Hamiltonian, and Hamilton equations of motion. This 
formulation is equivalent to the saddle point discussion 
in the Martin-Siggia-Rose approach in paper II. 

(ii) The ensuing canonical phase space formulation al- 
lows for a discussion of the time-dependent probability 
distributions, i.e., the weak noise solutions of the Fokker- 
Planck equation, in terms of phase space orbits on con- 
served energy surfaces, governed by Hamilton equations 
of motion. The action associated with an orbit plays the 
role of a weight function in much the same way as the 
Hamiltonian entering the Boltzmann factor in the de- 
scription of thermodynamic equilibrium. In the kinetic 
nonequilibrium problem defined by the Langevin equa- 
tion the dynamic action yields the probability distribu- 
tions. 

(iii) In the canonical phase space formulation the un- 
derlying stochastic nature of the Langevin equation and 
the relaxational character of the solutions of the Fokker- 
Planck equation are reflected in the topological struc- 
ture of the energy surfaces. A structure which differs 
markedly from the energy surface topology for ordinary 
dynamical problems In particular, the zero-energy man- 
ifold which determines the stationary state, i.e., the sta- 
tionary probability distribution, has a two-fold subman- 
ifold structure, including a hyperbolic stationary point, 
which in the simple case of a single stochastic variable 
corresponds to the unstable maximum of an inverted po- 
tential. Moreover, the waiting time for the orbits pass- 
ing close to the stationary point accounts for the Marko- 



3 



vian behavior of the probability distributions. Finally, 
the long-time orbit close to the zero-energy manifold de- 
termines via the action the time-dependent probability 
distribution. 

(iv) In the case of a few degrees of freedom the canon- 
ical phase space approach yields the established results 
following from an analysis of the Fokkcr-Planck equa- 
tions. On the other hand, in the case of many degrees 
of freedom, i.e., the field theoretical case, the Fokker- 
Planck equation becomes an unwieldy multi- dimensional 
differential-integro equation, and the canonical phase 
space approach, replacing the Fokkcr-Planck equation (in 
the weak noise limit) with coupled canonical field equa- 
tions yields, in addition to providing an alternative point 
of view of the stochastic processes in terms of dynamical 
system theory, a methodological advantage; particularly 
in the case where we can determine the zero-energy man- 
ifold explicitly. 

(v) In the field theoretical cases of the noisy Burgers 
or KPZ equations in one dimension, we can for special 
reasons identify the zero-energy manifold and determine: 
(1) the stationary distribution, (2) the long-time diffu- 
sive mode contribution to the time-dependent skew dis- 
tribution and (3) a heuristic expression for the short-time 
(transient) soliton mode contribution. In the interesting 
case of higher dimensions we are only able to make some 
general statements. 

The paper is organized in the following manner. In Sec 
II we consider the generic case of a nonlinear Langevin 
equation for many stochastic variables driven by white 
noise. We analyze the associated Fokker-Planck equation 
in the weak noise limit and set up the canonical phase 
space formulation. In Sec. Ill we consider as an example 
the stochastic overdamped motion in a harmonic poten- 
tial. In Sec. IV we apply the formulation to the Burgers 
and KPZ equations and derive expressions for the distri- 
butions. Finally, in Sec V we present a discussion and a 
conclusion. 



II. THE CANONICAL PHASE SPACE 
APPROACH 

Path integral formulations of the Fokker-Planck equa- 
tion in the field theoretical case and aspects of the canon- 
ical structure have been discussed in the literature, see 
|f37f for further referencing, we believe, however, that the 
present emphasis on the canonical phase space formu- 
lation as a practical tool is new. For this purpose we 
here set up the general canonical phase space formalism. 
Adhering to the notation in |p7fl we consider a general 
Langevin equation with additive noise [p8|,p9[ , 



-F n (q„ 



Vn 



(2-1) 



dq n _ _1 
dt ~ 2' 

Here q n , n = 1, ..N, is a set of time-dependent stochastic 
variables. The index n is discrete but is readily general- 
ized later to the field theoretical case of infinitely many 



degrees of freedom where n typically includes the spa- 
tial variables. The forces F n (q m ) are general functions of 
q n . In the linear case of coupled (overdamped) oscilla- 
tors, F n — 2Q nm q m (t), where f2 nm is a damping matrix. 
Finally, the equation is driven stochastically by a white 
noise term rj n with a Gaussian distribution and correlated 
according to 



(Vn{t)Vm(t')) = AK nm 6(t-t') 



(2.2) 



Here K nm is a constant, symmetric, positive-definite 
noise matrix of 0(1) and the correlations are character- 
ized by the noise strength A. 

Introducing the notation V„ = d/dq n the Fokker- 
Planck equation for the (conditional) prob ability distri- 
bution P(q n , t, q' n ) associated with Eq. ( |2.l| ) has the form 
"31 



dP 1 , 

-qT = ^4F n P + AK nm V m P] 



(2.3) 



including a drift term V n (F n P), arising from the deter- 
ministic force F n , and a diffusion term Ai^„ m V„V m P, 
originating from the noise rj n . 

In the equilibrium case, choosing K nm = S nm and set- 
ting F n = V„$, corresponding to an effective fluctuation- 
dissipation theorem and an underlying thermodynamic 
free energy $, Eq. (2^3) admits the stationary solution 
P s t oc exp [— $/A] with A entering as a temperature and 
mathematically as a singular parameter, i.e., the Boltz- 
mann distribution. Using this as a guiding principle we 
searc h in the general nonequilibrium case for solutions to 
Eq. (U) of the form 



P oc exp 



A 



S 



(2.4) 



where the weight function S replaces the free energy $ 
in the equilibrium case. By insertion and keeping only 
terms to leading order in A, it is easy to show that S sat- 
isfies an equation of the Hamilton- Jacobi form KG K32] , 



— + H(q n ,V n S) = 0, 
at 



(2.5) 



where, introducing the canonical momentum and energy 

Pn= V„S , (2.6) 
E= H , (2.7) 

the conserved energy or Hamiltonian is given by 

H = -K nm p n p m - -F n p n . (2.8) 

From the symplectic structure and dynamical system the- 
ory the canonical phase space structure follows immedi- 
ately. The action S has the form ]6C|] 



S 



(2.9) 
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and from the ensuing principle of least action we derive 
the Hamiltonian equations of motion, dq n /dt = dH/dp n 
and dp n /dt — —dH/dq n , 



dq n _ K 
at 



dt ~ 2 Pm n m ' 



(2.10) 
(2.11) 



for the orbits in p n q n phase space. 

The above formulation allows a simple interpretation 
of the solution of the Fokker- Planck equation ( [2.3| ) in the 
weak noise limit A — ► in terms of orbits in a canonical 
phase space. In order to determine the transition prob- 
ability P(q n ,T,q' n ) for a configuration q' n at t — to a 
configurati on at t — T, w e simply solve the Hamilton 
equations ( |2.1C| ) and ( 2.11 ) for an orbit from q' n to q n 
traversed in time T an d, subsequently, evaluate the ac- 
tion according to Eq. (2.£), yielding the weight function 
in Eq. ([Oh, i.e., 



P(q n ,T,q n ) oc exp 



1 

A 



T,q n 



o,q' n 



(2.12) 



We notice that the relevant orbit is determined by the 
initial and final values q' n and q n and the elapsed time 
T. The canonically conjugate mom entum p n is a slaved 
variable determi ned by Eq. ( |2.1l| ) and parametrically 
coupled to Eq. ( 2.10 ). Also, unlike the case in ordinary 
mechanics, the energy E = H in Eq. ( |2.S| ) is not the cen- 
tral quantity in the present interpretation. The traversal 
time T is the important variable and the energy manifold 
E(T) on which the orbit from q' n to q n lies is a function 
oiT. 



The stochastic nature of the Langevin equation (2.1) 
and the properties of t he weak noise solution of the 
Fokker-Planck equation (2^) are reflected in the topolog- 
ical submanifold structure of the energy surfaces in p n q n 
phase space. Unlike an ordinary mechanical problem H 
is not bounded from below and does not separate in a 
kinetic ener gy a nd a potential energy only depending on 
q n . In Eq. (2.8) the potential —(l/2)F nPn is momentum 
(velocity) dependent and gives rise to unbounded motion. 
Assuming for simplicity that F n — > for q n — > the en- 
ergy surfaces have the submanifold structure depicted in 
Fig 1. 

The origo in phase space constitutes a hyperbolic sta- 
tionary point, that is a saddle point determined by the 
zero-energy submanifold p n = and the zero-energy sub- 
manifold defined by K nm p m — F n being orthogonal to p n , 



K 



nmPin 



F n ± Pn 



Assu ming F n (q m ) ~ 2f2„ m <7 m for small q n the Hamil- 



tonian (2.8) is quadratic in p n and q n and a stability 
analysis can easily be carried out. In accordance with the 
present physical interpretation we assume that the sta- 
bility or damping matrix f2 nm implies an unstable p n = 
submanifold and a stable submanifold K mn p m — F n -L p n - 



The orbits in phase space close to the zero-energy man- 
ifold are thus those depicted in Fig. 1. In the harmonic 
oscillator picture which applies close to the stationary 
point this behavior corresponds to the motion in an in- 
verted parabolic potential as discussed in more detail in 
Sec. III. 




FIG. 1. Canonical phase space in the general case. The 
solid curves indicate the zero-energy transient submanifold (I) 
and stationary submanifold (II). The stationary saddle point 
is at the origin. The finite time (T) orbit from q' n to q n mi- 
grates to the zero-energy submanifold for T — > oo. 

The stationary state is given by orbits on the zero- 
energy manifold whose structure thus determines the na- 
ture of the stochastic problem. Assuming that E(T) oc 
exp [— const. T] for T — > oo the stationary state 



P s t(q n ) = lim P(q, 

T — >oo 



. T » q'n) 



is obtained from Eq. ( |2 . 1 2| ) , i.e., 

Pst(q n ) oc exp 



1 r , dq n 

A Ir, dtPn ^T 



(2.13) 



(2.14) 



We note that for T — > oo the orbit from q' n to q n con- 
verges to the zero-energy manifold, i.e., P s t is determined 
by the infinite-time orbits on the zero-energy manifold. 

The basic structure of phase space, depicted in Fig. 1, 
allows for a simple dynamical discussion in terms of dy- 
namical system theory Jl8| , |62| , |63[ of the approach to the 
stationary state of a damped noise-driven stochastic sys- 
tem. We first consider an orbit on an E = surface from 
q' n to q n in time T. The energy surface E(T) depends 
on T and in the limit T — > oo, E — + in order to attain 
the stationary state. For E — > the initial part of the 
orbit moves close to the p n = submanifold and from 
Eq. (2.10) is determined by 



dqn 
dt 



1 



-Fn 



(2.15) 



i.e., the deterministic noiseless version of the Langevin 
equation fl2.1| ). In the absence of noise the motion is 
transient and damped. Near the transient submanifold 
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Pn = the corresponding action S ~ and the proba- 
bility P = const., corresponding to a deterministic be- 
havior. The orbit slows down near the stationary point 
in phase space before it picks up again and moves close 
to the other stationary submanifold K nm p m — F n _!_ p n . 
This part of the orbit carries a finite action S 1 , i.e., P de- 
pends on q n , terminates in q n at time T, and corresponds 
for T — > oo to the stationary state. The Markovian be- 
havior, i.e., the loss of memory or the independence of the 
initial configuration q' nl is thus associated with the long 
(infinite) waiting time near (at) the stationary point. 

Whereas the transient submanifold p n — 0, yield- 
ing E = 0, is consistent with the Hamiltonian equa- 



tions (2.1C) a nd ( 2T1) and gives rise to the determin- 
istic equation ( 2.15 ), the other possibility of imposing a 
zero-energy submanifold by setting K n m P m = Fn will 
in general violate Eqs. ( 2.10 ) and ( 2.11 ). Assuming for 
simplicity K nm = 5 nm (note that the symmetric noise 
matrix K nm can always be diagonalized by a suitable 
choic e ofp w) we obtain, inserting p n = F n in Eqs. ( 2.10| ) 
and (2.11), that the relationship F m (V m F n — V n F m ) = 
must hold. In the special case where F n = V ra $, corre- 
sponding to an effective fluctuation-dissipation theorem 
and an underlying free energy the above constraint is 
trivially satisfied and we obtain the time-reversed equa- 
tion of motion dq n /dt = (l/2)F n governing the orbit on 
the p n — F n stationary submanifold. It is interesting to 
notice that the damped motion on the transient p n = 
submanifold is precisely equal to the time-reversed grow- 
ing motion on the stationary p n — F n submanifold. Fi- 
nally, it is an easy task to determine the stationary dis- 
tribution 



P st oc exp 



A 



(2.16) 



by insertion of p n — V n $ in Eq. (2.14) and integrat- 
ing over time., in agree men t with the solution of the 
Fokker-Planck equation ( |2.3| ) in the stationary case set- 
ting K nm = 5 nm and F n = V„$. 

We are led to the conjecture that for systems in ther- 
mal equilibrium where the force F n is derived from a 
thermodynamic free energy $, F n = V„$, the infinite- 
time orbits on the stationary zero-energy submanifold 
actually converge to the sub mani fold p n — F n , yielding 
the stationary distribution ( [2.16 ). In the general case 
of a driven stochastic system with a force F n not deriv- 
able from a free energy the only constraint is given by 
KnmPm — F n _L p n and we ha ve to solve the Hamil- 
ton equations ( 2.10| ) and (2.11) in order to determine 
the phase space orbits. 

We finally wish to comment on the connection between 
the present canonical phase space approach and the for- 
mulation in terms of a Martin-Siggia-Rose (MSR) path 
integral presented in paper II. For simplicity we consider 
only the case of a single stochastic degree of freedom q 
and set the noise matrix K nm = 1. 

All the relevant properties are extracted from the gen- 
erator 157P 



Z(/i) = (exp[t / dt fi(t)q(t)}) , 



(2.17) 



where (• • •) denotes an average with respect to the Gaus- 
sian noise distribution 



P{rf) oc exp 



1 

2A 



dt r]{tf 



(2.18) 



The Langevin equation (2.1) for one degree of free- 
dom is enforced by the delta function constraint 
/ JJ t dt,J5(dq/dt + (l/2)F — rj) = 1, where the Jacobian 
J = exp [(1/4) J dtdF/dq] |57||. Exponentiating the con- 
straint and in the process introducing an additional noise 
varia ble p, averaging over the noise n according to Eq. 
( |2.18| ), and scaling p, p p/A, we obtain the expression 



Z(fi) cx J Y\_dpdqexp 



exp 



i J dt n{t)q{t) , 
(2.19) 

where the action has the Feynman form pj,pq| 



•^'msi; : / <il ( p-Tj; — HmSR 



with complex Hamiltonian 



H 



MSI! 



1 t? i 2 

= — pb p 



.AdF 
'J~dq 



(2.20) 



(2.21) 



For the probability distribution P(q, t, q') we find in par- 
ticular 



P(q, t, q') oc / Y\ dpdq exp 

Jq> t 



(2.22) 



where 5*msr = J dt{pdq/dt — Hmsr) and the path inte- 
gral samples all orbits from q' = q(Q) to q = q(t) weighted 
with S'mrs; note that the noise field ranges freely. 

Reconstructing the underlyin g qu antum mechan- 
ics, yielding the path integral ( |2.22 ), P(q,t,q') can 
be regarded as a matrix element of the evolution 
operator exp [— iHMSRt] m a q-basis, P(q,t,q') oc 
(q\ exp (— iHMSB,t)\q') , where -ffMSR is the quantum ver- 
sion of Hmsr.- It follows that P(q, t, q') then satisfies the 
Schrddinger equation 



dP 



(2.23) 



with A playing the role of an effective Planck constant. 
The noise variable p becomes the momentum operator 
p = —iAd/dq and we obtain, inserting p — > p in Eq. 



H 



1 



MSR 



dq 2 



i(pF) ordcl + 



AdF 
~2~dq 



(2.24) 



6 



where since [p, q] = —iA and F depends on q we still 
have to specify the o rderin g in (pfLdcr- 

Comparing Eqs. fl2.23p and (|2.24j ) with the Fokker- 
Planck equation (2.3) in the present case, 



dP 
~dt 



dq 2 



F- 



dq 



dF 
dq 



P 



(2.25) 



we find agreement provided we choose the symmetric or- 
dering (pi 7 ') or dcr = \ {pF+Fp) . Alternatively, we are free 
to choose a normal ordering (pF) OT( ] LCr — pF and n eglec t 
the Jacobian contribution (A/2)dF / dq in Eq. ( [2.24 ). 
The Fokker-Planck equation then becomes the under- 
lying Schrddinger equation for the path integral with a 
complex non Hermitian Hamiltonian, 



Hfp = i — 



A 2 — + A—F 

dq 2 dq 



(2.26) 



The non Hermitian form of H pp with the p operator on 
the left ensures that Eq. (2.25) has the form of a conser- 
vation law ensuring the conservation of probability. 

In the limiting case F — 2uq the Hamiltonian Hmsr 
describes a harmonic oscillator and it is easy to see that 
the Jacobian contribution i(A/4)dF/dq — i(A/2)u in 
-Hmsr precisely cancels the zero point motion and en- 
sures a stationary state for t — > oo. 

Finally, in the cl assic al weak noise limit for A — > 0, 
the path integral (2.22) is dominated by the classical 
or stationary orbits following from the principle of least 
action 8Sm.sk — and determined as solutions of the 
Hamilton equations of motion: dq/ dt = dHyisnJ dp and 
dp/dt = — OHusk/ d<l] the distribution being given by 
P oc exp [(«/A)5msr]- This procedure is, however, en- 
tirely equivalent to performing the limit A — > directly 
in the Fokker-Planck equation, yielding the Hamilton- 
Jacobi equation (2.5) and the present canonical phase 
space approach. 



III. THE HARMONIC OSCILLATOR 
EXAMPLE 



AN 



In order to illustrate the phase space method we here 
apply it to the simple case of an overdamped harmonic 
oscillator described by the Langevin equation and noise 
correlations 



dq 

dt = - UJq 



n 



(r,(t)r,(t')) = A6(t-t') , 
with associated Fokker-Planck equation 



dP 
~dt 



1B_ 
2dn~ 



A 



dP 
dq 



2uP 



(3.1) 
(3.2) 

(3.3) 



This system is well-known and easily analyzed p7| . |58[ 
The time-dependent pro bability distribution P{q, T, q') 
the solution of Eq. (|3.3|), is given by 



P(q, T, q) cx exp 



u> [q — q' exp {—bjT)Y 



A 1 - exp (-2uT) 
approaching the stationary distribution 

1 



(3.4) 



P s t(q) oc exp 



--u,q 



(3.5) 



in the limit T — » oo. 

We now proceed to derive these results within the 
canonical phase space formulation. Since F = 2toq and 
K = 1 for o ne degree of freedom, we obtain from Eqs. 
(^8|) - (|2.1l|) the Hamiltonian 



TT 1 2 



the action 



H 



(3.6) 



(3.7) 



and the Hamilton equations of motion 



dq 



Tt = p -" q 

dp 



(3.8) 
(3.9) 



The phase space is depicted in Fig. 2 and corresponds 
to the vicinity of the stationary point in Fig. 1 for one 
degree of freedom. 

For a single degree of freedom we can explicitly de- 
termine both zero-energy submanifolds: p — and 
p = 2ujq, and determine the orbits. On the transient 
p = submanifold dq/dt = —Luq, i.e., the determinis- 
tic equation of motion for rj = 0, with a damped solu- 
tion q — qo exp [— uxt], qo — q(0), corresponding to the 
damped orbit approaching the stationary saddle point at 
(p>S) = (0,0). The action associated with this orbit is 
S = 0, i.e., P — const., characterizing the determinis- 
tic motion. On the stationary submanifold p = 2ujq the 
Hamiltonian equations coincide, dq/dt = 2ujq, and we 
obtain a growing solution q = qo exp [cot], qo = q(0), asso- 
ciated with the orbit emerging from the stationary point. 
As discussed in Sec. II we note that the stationary orbit 
is the time-reversed mirror of the transient orbit. 

The complete solution of Eqs. (3.8) and (3.9) is also 



easily obtained. For an orbit from q' to q" in time T and, 
noting that p is a slaved variable, we obtain 



*(*)= 



q" sinh ojt + q' sinh w (T - t) 



p(t)= UJ 



sinhwT 
q"e wt - q'e-^-^ 
sinh luT 



(3.10) 
(3.11) 
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FIG. 2. Canonical phase space plot in the case of an over- 
damped oscillator. The solid curves indicate the zero-energy 
transient submanifold (I) and stationary submanifold (II). 
The stationary saddle point is at the origin. The finite time 
(T) orbit from q' n to q n migrates to the zero-energy submani- 
fold for T — *• oo. We have also indicted the sign of the energy 
E in the four domains. 

For large T the noise variable p is initially close to 
zero, corresponding to the transient deterministic regime; 
for t close to T p eventually leaves zero and approaches 
the limiting value 2uq" , corresponding to the stationary 
regime. Likewise, we note that q ~ for T — > oo for 
most t. The behavior of q and p is depicted in Fig. 3. 



stationary point is then defined by the limiting orbit for 
T — > oo. 

Moreover, from the Hamilton equations and fl3.9| ) 
we readily deduce the second order Newton equation of 
motion 



fq 
alt 2 



■ u> 2 q 



(3.13) 



describing the orbit in an inverted harmonic potential 
— (l/2)oj 2 q 2 , and allowing for a simple discussion of the 
motion in pq phase space. 

The finite energy orbits fall in two categories depend- 
ing on the sign of E. For E > the orbits pass through 
the unstable maximum of the inverted potential with fi- 
nite momentum; for E < the unbounded orbits are 
confined by the potential to either positive or negative 
values of q. The limiting case E = corresponds to an 
orbit approaching the maximum, the hyperbolic station- 
ary point, with zero momentum. This point represents 
an unstable equilibrium where the particle spends an infi- 
nite amount of time, corresponding to the establishment 
of Markovian behavior, i.e., the loss of memory and in- 
dependence of the initial configuration q" . The motion 
is depicted in Fig. 4. 



> 2 q 2 



(a) 



♦ t 



(b) 



q' q" 



2coq" 



FIG. 3. Plot of q and p as functions of t in the case of 
an overdamped oscillator. In a) we depict the dependence of 
q; for large T the coordinate q stays close to the stationary 
saddle point. In b) we show the dependence of p; for large T 
the momentum p is initially close to the transient submanifold 
p = but eventually moves on to the stationary submanifold 
p = 2uq. 

The orbit from q' to q" traversed in time T lies on the 
energy surface given by 



E = 



u> 2 q" 2 



q ' 



2q'q" coshwT 



sinli loT 



(3.12) 



For fixed q' and q" the energy E is a function of T, 
E = E(T). In the limit T -> oo, E -> and the or- 
bit converges to the zero-energy manifold as indicated in 
Fig. 2. The actual T = oo orbit passing through the 



FIG. 4. In the case of the overdamped oscillator the orbits 
in pq phase space depicted in Fig. 2 corresponds in q space to 
the motion in an inverted parabolic potential. The unstable 
maximum corresponds to the stationary saddle point. 



In terms of the explicit solution ( 3.10 ) and ( 3.11 ) we 
finally derive the action associated with the orbit, 



o = (jJ- 

and recover from P 
and in the limit T - 



(q" - q'e~ ul ) 



T\2 



1 



-2uT 



(3.14) 



oc exp [-S/ A] the time-depen den t 
oo, stationary distributions (3.4) 



and (3.5), respectively. 

Th e Hamiltonian ([T^) and the equations of motion 
(|3.8| ) and (3.E), yielding the canonical phase space rep- 
resentation of the Langevin equation for a noise-driven 
overdamped harmonic oscillator, have the same struc- 
ture as the dynamical description of an ordinary har- 
monic oscillator. Shifting the momentum, p — » p + ujq, 
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H — > (l/2)p 2 — (l/2)wq 2 , describing the motion in an 
inverted harmonic potential as discussed above. The 
equations of motion now take the form dq/dt = p and 
dp/dt = Luq with solutions given as linear combina- 
tions of a g rowi ng and a damped solution as in Eqs. 
( |3 .10 ) and ( 3.11 ). However, performing a Wick rota- 
tion t — ► it = t in combination with the transformation 
to a complex momentum p — > — ip, note that p is basi- 
cally a dummy variable representing the noise, we obtain 
H — > —H osc , where i/ osc is the oscillator Hamiltonian 
#osc = 5P 2 + hu 2 q 2 , yielding the equations of motion 
dq/dr — p, and dp/dr = —cuq with bounded periodic 
motion in imaginary time r. The energy E = H osc is 
positive and the finite energy orbits in pq phase space 
move on concentric ellipses. The zero-energy manifold 
corresponds to the origin (p, q) — (0, 0) in phase space. 
In fact, subject to the Wick rotation t —> it and a rotation 
of p, p — > — ip, this phase space structure is mapped to 
the phase space structure in Fig. 2 with a hyperbolic sta- 
tionary point and unbounded orbits. As pointed out in 
Sec. II it is precisely in the energy surface topology that 
the stochastic problem differs from an ordinary dynami- 
cal problem; here exemplified in the context of the oscil- 
lator. Correspondingly, the action ([O]) for an orbit from 
q' to q" in time %T transforms to —iS osc + (w/2) (q" 2 - a' 2 ) 
where S osc is the action for a harmonic oscillator |&i| , |65f| 



5V>sr. 



2 sinajT 



[{q" 2 + q' 2 )cosuT-2q'q"] , (3.15) 



yielding P(q" , iT, q') in accordance with with Eq. (p.4| ). 

Finally, we present a simple calculation of the lead- 
ing correction at long times to the stationary distribu- 
tion (3J3) which will prove useful in our discussion of 
the Burgers-KPZ equations in the next section. For 
T — > 00 the orbit in phase space is close to the sta- 
tionary zero-energy submanifold p = 2ujq. Replacing 
the orbit from q' to q" at long times T with an orbit 
on the stationary manifold we obtain a constraint which 
allows us to simply evaluate the correction to P s t- Con- 
sequently, inserting the zer o-en ergy constraint p = 2cuq 
in the canonical equation ( |3.§| ) yields dq/dt — ujq with 
solution q" = exp (tvT)q' for q' and q" on the p = 2ujq 
manifold. Inserting this solution in the action (3.7) we 
obtain 



S = u{q" 2 -q' 2 ) 



ujq" z (l - exp (-2wT)) , (3.16) 



in accordance with an expansion of the exact result (3.14) 
to leading order in exp (— u>T). 



A. The general case 

First applying the canonical formulation to the Burg- 
ers equation (1.15) the index n in Sec. II now comprises 
both the continuous spatial coordinate x n , n = 1, ..d and 
the vector index of the slope field «„, n = 1, ..d, i.e., n — * 
x n ,n. Furthermore, we choose the noise matrix K nm 
and the forces F n according to the prescription, K nm — » 



-2(vV 2 u n + Xu m \7 m u 



V 2 n„ S(x n - x' n ) and F n 
With the identification q n u n (x m ) and p„ — > p n (x m ) 
we thus obtain the Burgers action and Hamiltonian den- 
sity 



■% C d d xdtf Pn ^-n B 



and the ensuing Hamilton equations of motion 

d_ 
di 

d_ 

at 



(4.1) 



W B = Pn [ U n + Xu m V m U n - -Vn^mPm ) , (4.2) 



Au m V m J u n = vV~u n - V„V m p ri 
Au m V m I p n = -vV 2 p n 



(4.3) 



Um) . (4.4) 

The time-dependent probability distribution is then in 
the weak noise limit given by 



P(u n ,T,u' n ) oc exp 



SB(u n ,T, u' n 



(4.5) 



The Hamilton equations (4.3) and (4.4) determining the 
orbits in p „n ra phase space thus replace the noisy Burgers 
equat ion ( |l,15| ) in the weak noise limit and the distribu- 
tion ( |4.5| ), evaluated for an appropriate orbit from u' n to 
u n traversed in time T, constitutes a weak noise solu- 
tion of the Fokker-Planck equation associated with the 
Burgers equation, 

dP(u n ,t) = 

at 

/s 
d d x- — [{vV 2 u n + Xu 7n VmU n )P{u n , t)] 
ou n 

V 2 Y[5{x n -x' x )P(u n ,t) 



+ 2 ./ ^ X Su n SuL 



(4.6) 



IV. CANONICAL FORMULATION OF THE 
BURGERS - KPZ EQUATIONS 

In this section we apply the canonical phase space 
method de velop ed in Sec. II to the Burgers and KPZ 
equations ( 1.15| ) and ( 1.17 ). 



The time-dependent and stationary distributions are de- 
termined by the orbits near and on the zero-energy mani- 
fold. From the general discussion in Sec. II it follows that 
the zero-energy manifold has a submanifold structure 
with a transient p n — submanifold, a hyperbolic sta- 
tionary point at (u n ,p n ) = (0,0), and a stationary sub- 
manifold defined by i/V 2 u„ + Aw m V m u„ — (l/2)V„V m p m 







orthogonal to p n ; here treating the integration over x in 
J d d xH.B as an in ner product. On the transient subman- 
ifold p n = Eq. (i.4) is trivially satisfied and the orbits 
are governed by the noiseless Burgers equation 



9 x 
at 



(4.7) 



which is ana lyzed by means of the Cole-Hopf transfor- 
mation (1.2) Un — V„/i, h = (2^/A)lnw, with w sat- 
isfying Eq. (1.3) and solved by means of the Green's 



function (1.5) generalized to the <i-dimensional case. On 
the other hand, on the stationary submanifold defined 
by / d d xli.B = 0, determining the stationary distribution 
Pst (« n ) , t he o rbits are given by the coupled equations 
( |4.3| ) and (4.4) and will be discussed in the next section. 

For later reference we also pr esent the canonical for- 
mulation of the KPZ equation ( 1.17| ). Here we choose 

q n -» h(x), p„ -> p(x), K nm ~> Iln^C^n _ X 'n)> and 

F n — > — 2(i^V 2 /i + (A/2)V„W„/i), and obtain action, 
Hamiltonian density, and equations of motion, 



Skpz= y d rf :rdt ( p— - 7i K pz 



A 



1 



^KPZ= P UV'll + -V n hV n h + - P 

dh A_ ,_ , 

at 2 



<9p 

at 



= -i^V 2 p + AV„(pV„/i) 



yielding the weak noise distribution 



P(/i, T, /i ) oc exp 



S K pz(h,T, ti) 



(4.8) 

(4.9) 

(4.10) 
(4.11) 

(4.12) 



as solution of the Fokker - Planck equation 
dP(h,t) 



dt 



dx Tk 



(vV 2 h+^V n hV n h)P(h,t) 



2 / ShSh' 



l[S(x n -x' n )P(h,t) 



(4.13) 



The KPZ formulation is, however, comp lete ly equiva- 
lent to the Burgers formulation. In Eq. (4.3) only the 
longitudinal component of the noise field p n couples to 
the longitudinal slope field u n = V„/i and we can with- 
out loss of generality as sum e that p n = V n is purely 
longitudinal since Eq. (4.4) is linear in p„; this prop- 
erty reflec ts the conserved noise V n ? 7 dri ving t he B urgers 
equation (1.15). Comparing Eqs (4.3) and (4.4) with 
Eqs. (4.10) and (4.11) we obtain complete equivalence 
by choosing V m p m = —p or p = -V 2 0. 



B. The one dimensional case 

In one dimension and focussing on the slope field u, 
which in many respects is the natural variable in dis- 
cussing a growing interface, the canonical equations (4.3) 
and (4.4) take the simple form 



- AuVj u= vV 2 u - V 2 p , 

- AmV ] p= ~vV 2 p , 



dt 
d_ 

dt 

originating from the Hamiltonian density 
Hb=P (vV 2 u + XuVu - -V 2 p 



(4.14) 
(4.15) 

(4.16) 



We note that both u and p are sc alar fields and that the 
A-dependent term on the RHS of (4.4) cancels. Also, sub- 



ject to th e shift transformation p = v(u — ip) Eqs. ( 4.14 ) 
and ( 4.15| ) are identical to the equations (1.10) and (1.11) 
discussed in paper II. 

It is an important property of the one dimensional case 
that we can determine the explicit form of the station- 
ary zero-energy manifold, as was the trivial case for one 
degree of freedom discussed in Sec. II. For 



p = 2vu 



(4.17) 



the canonical equations ( t4.14[) a nd (4.15) become identi- 
cal and the energy density ( |4.16j ) takes the form of a total 
derivative, Hb —> (2/3)A^Vu 3 , yielding a vanishing total 
energy Eb — J Hb = for vanishing slope field at the 
boundaries. Owing to the vector character of u n and p n 
and the presence of the X{p n V m u m — p m X7 n u m ) term in 



Eq. (4.4) such a transformation does not seem possible 
for d > 1 and we dont have a similar contraction of the 
stationary submanifold. 

In other words, in d = 1 the orbit from u' to u" in 
time T for T — > oo does not only approach the zero- 
energy submanifolds p = and vV 2 u + XuVu — ^V 2 p 
orthogonal to p but actually converges to the submani- 
fold p = Ivu on the stationary submanifold. This phase 
space behavior is depicted in Fig. 5. In Fig. 5a we show 
the contraction of the stationary manifold. In Fig. 5b we 
depict the orbits in pu phase space in a similar manner 
as in Fig. 2. 

We finally wish to present a plausibility argument for 
the a ttrac tion of the orbits to the submanifold given by 
Eq. (4.17). Denoting the deviation from the su bman ifold 
by 8u and inserting p = 2v(u + Su) in Eqs. (4.14) and 
(4.15) we obtain to leading order in Su, 



0_ 

Ot 



AwV 6 



^ 2 5u . 



(4.18) 



Noting that d/dt — AitV in invariant under the Galilean 
transformation x — > x — Xugt, u — > u + uq and choos- 
ing a local frame with vanishing it, the Fourier modes 
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Suk oc exp [— vk 2 t\ — > for large t, implying that the 
orbits approach the zero-energy submanifold. 



p=2nu 





FIG. 5. Here we depict the phase space behavior in the 
case of the noisy one dimensional Burgers equation. In a) we 
show the contraction to the zero-energy submanifold p — 2i/u, 
characteristic of the one dimensional case. In b) we show sim- 
ilar to Fig. 2 the orbits in pu phase space. 



V. DISCUSSION AND CONCLUSION 

In this final section we derive results for the probability 
distributions and attempt to draw some general conclu- 
sions on the basis of the canonical phase space approach 
to the noisy Burgers and KPZ equations presented in the 
previous sections. 



A. The one dimensional case 



The timc-dcpcndent distribution (4.5) is dete rmined 
by the form of the action Sb(u,T) in Eq. (hi). The 
following analysis implies that the action has the generic 
form 



Sb = S st (u) + S d is(u, T) + S so i{u, T) , 



(5.1) 



where S^t yields the stationary distribution (1.9), Sdia 
gives rise to corrections due to the linear diffusive modes, 
and S'sd originates from the nonlinear soliton modes; 
both Sdiff and S so \ must vanish in the limit T — > oo so 
that we attain the stationary distribution given by S s %. 

In the linear Edwards- Wilkinson case for A = only 
diffusive modes contribute and there is no growth. In 
wave number space the field equ ations (]4. 14 ) and ( 4. 15] ) 
take the same form as Eqs. (3.8) and ( p?9 ) in Sec. III. A 
straightforward generalization of Eqs. (3.1C) and (3.14) 
then yields the orbit, u k = U}~(T), u k 
u k = vk 2 , 



Mfe(O), and 



u k {t) 



u' k sinhw/ct + u' k sinhcjfc(T — t) 
sinhajfcT 



(5.2) 



S\in = V 



dk \uk — u' k exp(— uikT)Y 
2^ 1 - exp(-2w fc T) 



(5.3) 



We note that in the limit T — > oo the action Sn n (u, T) — * 
v f(dk/2ir)\uk\ 2 i n ac cordance with the stationary dis- 
tribution in Eq. fll.9| ) expanded in wave number space. 
Since at long times Uk — exp(tOkT)u' k we also find the 
correction 



'diff 



dk 
2^ 



u fc | 2 exp(-2w fc T) 



(5.4) 



From Eqs. (5.3) and (5.4) we observe the simple scaling 
property, u -> fiu, S iin -> /J, 2 S lin , and S diS -> fJ. 2 S dlS , 
i.e., scaling the slope or height field with a factor fi the 
action scales with fi 2 . Th is be havior is com patible with 
the equations of motion (4.14) and ( |4.15 ) provided we 
scale p — » /ip. 

As regards the time dependence of Su a and Sdiff we 
identify the crossover time 



t: 



diff 



1 

vk 2 



(5.5) 



depending on the wave number k. In the thermody- 
namic limit L — > oo the wave number k has a continu- 
ous range and the crossover time diverges in the infrared 
limit k — > 0. Consequently, we do not have a separation 
of time scales. 

Since the saturation width of an interface is a finite 
size effect time scale separation only occurs for a finite 
system. In the present linear case this is associated with 
the quantization of the wave number k ~ 1/L; note that 
the n = mode is related to the global conservation of 
slope , i.e., J dx u is a constant of motion, and we have 



jrdiff ^ ^_ 



(5.6) 



and action, here Uk = Uk(T), 



From the general discussion of a growing interface, T co oc 
L z , where z is the dynamic scaling exponent, and we 
readily obtain z = 2 in accordance with the diffu- 
sive mode contribution with dispersion ui^ = vk 2 . For 
T <C the diffusive modes contribute to the time- 

dependent distribution, whereas for T 3> T~[ we cross 
over to the stationary regime, Sdiff fa; T ) — > 0, and we 
approach the stationary distribution ([L^) determined by 
S st {u). 



1. The stationary distribution 

In the one dimensional case the stationary distribu- 
tion is known p4] and has the symmetric Gaussian form 
given by Eq. ( [Of ) . This follows directly from the station- 
ary Fokker-Planck equation ( |4.6| ) where the A-dependent 
term for a Gaussian distribution becomes a total deriva- 
tive and thus yields a vanishing contribution; an argu- 
ment which only holds in one dimension. The slope 
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field u(x) is thus uncorrelated beyond a finite correla- 
tion length £ which is zero for the Burgers equation and 
microscopic for lattice models falling in the same univer- 
sality class. The height h(x) = f dxu performs random 
walk yielding according to Eq. (1.12) the roughness ex- 
ponent C = 1/2. 

Within the present canonical phase space formulation 
the stationary distribution follows immediately from the 
structure p = 2vu of the stationary submanifold. As 
in the harmonic oscillator case in Sec. Ill, the diffusive 
modes imply that Eb — > for T — ► oo. Thus inserting 
p = 2vu in the action (4.1) in the one dimensional case 
and performing the time integration, we obtain in the 
limit T —> oo 



POO Qu f 

S s t(u) = J dxdt p— — v J dx u{x) 2 



(5.7) 



and by insertion in Eq. (4.5) the distribution fll.9|) . 

Whereas the effective fluctuation-dissipation theorem 
valid in one dimension implies that the stationary distri- 
bution is Gaussian and symmetric in the slope u and in 
the height field h, measured relative to the mean height 
(h), the time-dependent distribution, converging towards 
the stationary one, is expected to exhibit an asymmetric 
shape corresponding to the predominance of peaks in h 
in the growth direction. 



2. The long-time skew distribution 

The phase space approach also allows us to estimate 
the long-time corrections to the stationary distribution. 
Following the reasoning in Sec. Ill we replace for large 
T the orbit near the submanifold p = 2vu with an orbit 
on the submanifold. Inserting p 



2vu in Eq. fl4.14| ) the 
orbits on the stationary submanifold are governed by the 
noiseless Burgers equation with v replaced by - v 1 



du 
dt 



— AuV I u = —is\7 2 u . 



(5.8) 



This equation is readily solved by means of the Cole-Hopf 
transformation (1.2) with solution given by Eqs. ( |1.4| - [i~G| ) 
with -v substituted for v. For the action ( |4.l| ) we then 



obtain Sb = v J dx[i 
probability distribution 



and for the time-dependent 



P{u,T) oc P st (u)P s kew(u,T) , 



(5.9) 



where the sym met ric stationary distribution P s t(u) is 
given by Eq. ( |l.9|) and the time-dependent skew cor- 
rection by 



P s kew{u,T) = exp 



£ Jdx u'(xf 



(5.10) 



with u' = V/i' and u = Vh related by the Cole-Hopf 
transformation 



exp 



2v 



h'(x) 



= J dx' G(x — x 1 , T) exp 



2;/ 



h(x') 



(5.11) 



Note that since / dx'G = 1 the correction u' 2 vanishes in 
the limit T — > oo and P s kcw — * 1- 

In order to examine the skewness of the distribution it 
is convenient to eliminate the stationary component by 
forming the ratio 



P{u,T) 



P(-u,T) eXP [^ / dx ^ U +- U '- ) \ ' (5 ' 12) 
where according to the Cole-Hopf transformation ( |5.11 ) 



exp 



-—h' ± 
2v ± 



— J dx' G(x — x, T) exp 



^2v 



(5.13) 



Inserting the Green's function (1.5), G(x,T) — 
[A-kvT] -1 / 2 exp[— x 2 /4j/T], we consider first a few simple 
cases. 

For a constant slope u = uq, i.e., h = u$x — ho, 
we obtain, performing the Gaussian integration, h'± = 
±{uqx + ho) — XTuq/2. We note that the growth term 
AwVu as expected gives rise to a time-dependent term 
in h'±. This term, however, is compensated for by trans- 
forming to a co-moving frame as in the KPZ equation 
(1.17). The slope it', however, is independent of T and 
we obtain it' ± = ±uq, yielding 5b = and P = const., 
i.e., no dynamics. This is consistent with the fact that 
u = uq an d p — po trivially satisfy the field equations 
( 4.14 ) and ( 4.15|) yielding Eb = 0, and thus corresponds 



to a stationary state, as also discussed in paper II. 

Choosing a slope depending linearly on x, u = 2sqx, 
corresponding to a parabolic height profile, h = sox 2 +ho, 
we obtain it' ± = ±2sox/(l ± 2ATso), yielding the skew- 
ness ratio 



P(u,T) 
P(-u,T) 



exp 



32 v\ 



TL 3 



3 A (l-(2ATs ) 2 ) 2 



(5.14) 



where we have introduced the size L of the system. The 
expression ( |5.14 ) only holds for ATso < 1, the important 
aspect is, however, the dependence on the sign of sq, i.e., 
the slope of the slope or the bias of the height profile. 
For So > corresponding to a parabolic shape of h with a 
minimum, i.e., a downward peak, P{u,T) / P(—u,T) < 1, 
whereas for sq < 0, yielding an upward peak in h, we 
have P(u,T)/P(-u,T) > 1. This behavior implies that 
the distribution is skew at finite times and that the up- 
ward peaks in h statistically are more pronounced than 
the downward peaks, i.e., the distribution is biased and 
changes asymmetrically towards the symmetric station- 
ary distribution. 

We can also gain some insight in the inviscid limit 
v —> from a saddle point calculation along the lines 
of similar treatments of the noiseless Burgers equation 
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|32| . |33| . F rom Eq. (5.11) we obtain, inserting the Green's 
function 



skew 



PoPx , 



exph^' ± ] - J dx' [4^T]- 1 /2 exp[ _ (1/2i , )0± ] 5 

(5.15) 



where (f>±(x,x') = (x - x') 2 /2T ± Xh(x'). In the limit 
v — > the integral is dominated by the local min- 
ima given by d<j)±(x 1 x') / dx' — and the condition 
d 2 (j)±(x,x')/dx' 2 > 0. The solutions x'± are thus de- 
termined by the implicit equation X _|_ ■ 3? — 
together with the conditions ±(du/dx') x i =x ' ± > — 1/AT, 
and we obtain the ratio 



P(u,T) 
P(-u,T) 



exp 



x' + ) 2 - (x-x'J 2 



(AT) 5 



(5.16) 



which can be used in order to analyze the skewness of var- 
ious profiles u(x). Note that P(u,T)/P(-u,T) -> 1 in 
the limit T — -> oo, corresponding to vanishing skewness. 
In Fig. 6 we have depicted the saddle point construction. 




FIG. 6. Here we show the saddle point construction valid 
in the inviscid limit v — > 0. The saddle point condition 
x'_ — x = XTu(x'_) determines x'_ as the intersection be- 
tween the line with slope 1/AT and the slope profile u(x'). 
The intersection of the line with the x' axis determines x. 

More insight into the dynamics underlying t he sk ew 
long-time distribution is gained by expanding (5.11) in 
the nonlinear coupling. Choosing a compact notation, 
U x = u(x), h x = h(x), G XX '{T) = G{x - x',T), and 
5xx' = S(x — x'), we obtain to leading order in A 

u' x = J dx' G xx <{T)u x > 

~{^/ v ) J dx' G xx '(T)u 

x' hx' 

+(X/v) J dx'dx" Gxx>{T)Gxx»{T)u x 'h x » . (5.17) 

Correspondingly, P s kew factorizes in a component inde- 
pendent of A and a component depending on A, 



P a = exp 
P A = exp 



V -\dxdx' G xA 2T)u x u x 

2A r 

-— J dxdx'dx" F xx ,. x , x „(T)u x u x ,h x 



(5.18) 
(5.19) 



(5.20) 



Here the kernel F is given by 

Fxx' ,x'x" (P) 

G xx i{2T)5 x i x n - J dy G yx (T)G yx <(T)G yx »{T) . 

(5.21) 

For vanishing A only Pq contributes. In wave 
number space we thus obtain, introduci ng U k — 
J dx exp (ikx)u(x) and noting from Eq. (1.5) that 
Gk(t) = exp [— i/k 2 t], the distribution 

PK,T)cxexp --J g | Ufc | 2 [l-exp(-2 I ,fc 2 T)] 

(5.22) 



This result is completely equivalent to Eq. (3.16) for the 
damped oscillator discussed in Sec. III. For A = 0, cor- 
responding to the Edwards Wilkinson case, only gapless 
diffusive modes Uk <x exp [~vk 2 t] contribute to the time- 
dependent distribution; we also note that the distribution 
remains symmetric. 

To first order in A the diffusive modes interact a nd w e 
obtain a correction to P s kew given by Pa in Eq. ( 5.20 ). 
This contribution is odd in u and characterized by the 
kernel F. In other words, the scattering of the diffusive 
modes on one ano ther due to the term AuVm in the Burg- 
ers equation (5^8) yields a skew distribution in u. In the 
limit T — ► oo this term vanishes and we obtain the sym- 
metric stationary distribution. Correspondingly, in wave 
number space we have to order A 



P\{u k ,T) cx 



2A fdkdk' 

-r / ——tk.k'\J-)ukU-k-k'hk' 

A J Z7T Z7T 



(5.23) 

F k ,k'{T) = G k (2T) - Gk{T)Gk+k-{T)Gk-{T) , (5.24) 

showing the interaction between the various fc-modes (the 
cascade). 



3. The short-time skew distribution 

At shorter times, i.e., T < T c ^ lff , the approximation 
of replacing the orbit near the zero-energy submanifold 
with an orbit on the submanifold ceases to be valid and 



we have to consider the equations of motions (4.14) and 
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(4.15) in more detail in order to identify the contribution 
to S so \- 

Although the noiseless Burgers equation is exactly sol- 
uble by means o f the Cole- Hopf transformation, the equa- 
tions of motion (4.14) and (4.15) describing the noisy case 
are presumably not exactly integrable. They do, how- 
ever, admit special permanent profile or solitary wave so- 
lutions with superposed linear diffusive modes. It more- 
over follows from the path integral formulation in paper 
II that an arbitrary interface profile can be represented 
by a dilute gas of solitons, at least in the inviscid limit 
for small v, where we can neglect soliton overlap contri- 
butions. 



In the form given by Eqs. (1.10) and (1-11) these equa- 
tions were discussed in detail in paper II where we iden- 
tified the elementary excitations. The spectrum consists 
of right hand and left hand nonlinear soliton modes with 
superposed linear diffusive modes. The localized soliton 
modes have a finite energy and thus correspond to the 
nearby phase space orbits approaching the zero-energy 
manifold. At long times the soliton energy must go to 
zero and the remaining superposed diffusive modes de- 
termine -Pskcw as discussed above. Note that the non- 
linear soliton mode can be regarded as a bound state of 
diffusive modes; this follows from the stability analysis in 
paper I and is a consequence of Levinson's theorem. In 
paper II we performed a shift transformation of the noise 
variable, p — > v(u — tp), i n ord er to express the Hamilto- 
nian density TL^ in Eq. (4.16) in a canonical form with 
a harmonic component describing the linear case, yield- 
ing the field equations ( l.lOj ) and (1.11). In the present 
context we summarize the soliton dynamics in terms of 
p and u in accordance with the present interpretation of 
the transition to the stationary state. The right and left 
hand solitons then play a different role in the weight of 
the interface morphology. In the static limit the soliton 
modes have the form 



u(x) = ±ito tanh 



2v 



(x - x ) 



(5.25) 



with amplitude uq and position x . Using the Galilean 
invariance of the field equations (4.14) and (4.15), i.e., ob- 
serving that the operator d/dt — AuV is invariant under 
the transformation x — > x — Xut, tt — > u + it, propagating 
solitons with boundary values u — > u± for x — > ±L, L 
is the size of the syst em, a re obtained by boosting the 
static solitons in Eq. ( 5.25| ). Moreover, the propagation 
velocity v is given by u+ and U- according to the soliton 
condition 



u+ 



2v 

T 



(5.26) 



which thus determines the kinetics and matching condi- 
tions for a multi-soliton configuration describing a grow- 
ing interface. 

The right hand soliton corresponds to p = and is 
according to Eq. ( |4.14 ) a solution of the noiseless Burg- 
ers equation (|l . l|) . Dynamical attributes are a feature 



of the noisy case and this soliton thus carries vanish- 
ing energy, Eb = f dxli.B — 0, vanishing momentum, 
IIb = J dxuVp = 0, and vanishing action Sb — 0, and 
corresponds to an orbit on the transient zero-energy man- 
ifold. We note that a single right hand soliton or a multi- 
soliton solution cannot satisfy the boundary condition of 
vanishing slope. 

The left hand soliton moves on th e sub manifold p = 
2vu and satisfies acco rding to Eq. (4.14) the noiseless 
Burgers equation ( |l.l| ) with v replaced by —v. It carries 
energy, momentum, and action given by 



Eb= o^ A ( u + ~ 



/ 2 2 
V[Uj_ — U_ 

Sb= -t;v\T\u + 



(5.27) 
(5.28) 
(5.29) 



Since u+ < u_ for a right hand soliton Eb is negative 
and Hb has according to Eq. (5.26) the same sign as the 
velocity v. The action Sb is positive and Galilean invari- 
ant. We also note that although the soliton is confined 
to the submanifold p — 2vu the energy is nonvanishing. 
This is associated with the nonequal bounda ry va lues u + 
and U- and also follows directly from Eq. (4.16), where 
inserting p — 2vu we obtain a total derivative. Hence, 
Eb = J dxHB — (2i/A/3)(itj_ — where u+ and U- 
are the boundary values. 
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FIG. 7. We depict the slope soliton morphology for a 
growing interface in a system of size L. In accordance with 
a growth situation we have imposed periodic boundary con- 
ditions. The left hand solitons carrying energy, momentum, 
and action are indicated by solid lines. 

As discussed in paper II the morphology of a grow- 
ing interface is determined by matching a set of right 
hand and left hand solitons according to the soliton con- 
dition fl5.26p . The right and left hand solitons are ex- 
act solutions of the damped and undamped (-v) noise- 
less Burgers equations, respectively. Their stability is 
associated with the nonequality of the boundary values, 
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corresponding to a nonvanishing slope cur- 



rent at the boundaries. In a multi-soliton configuration 
with vanishing boundary conditions current thus flows 
between the solitons. The current is generated by the left 
hand solitons and dissipated by the right hand solitons 
- this is another view of the cascade driving the noisy 
Burgers equation. The morphology is depicted in Fig. 7. 

The probability of a soliton morphology is determined 
by S so \(u,T). Assuming that S so \ — \F(u,T) a scal- 
ing argument similar to linear case, u — * fiu, p — > /ip, 
A — » fi X, and S so \ — » [i 2 S so \, foll owing from the gen- 
eral form of the equations of motion (4.14) and (4.15) and 
the action (4.1) impli es th at F — ► fi 3 F. This is consistent 
with the expression Q5.29 ) and we obtain 



s SO }(u,t) = g^Ary^ i 

lhs 



(5.30) 



where the summation is only over contributing left hand 
solitons (lhs). 

Beyond this point our discussion becomes necessarily 
more qualitative and heuristic since we dont possess a 
complete solution of the coupled field equations. The 
non-integrabili ty an d the constraint imposed by the soli- 
ton condition ( 5.26| ) imply that we only have available a 
dilute gas of right hand and left hand solitons. First we 
notice that fo r an infinite system the soliton velocity v 
given by Eq. ( |5.26 ) ranges freely, implying that the cen- 
ter of mass of the soliton amplitude u cm = (w+ + U-)/2 
also can take arbitrary values. Since T enters as a 
prefactor in Eq. ( |5.30 ) iS so i grows for T — > oo, i.e., 
P oc exp[— tSgoi/A] — > 0, and we are unable to identify 
a solitonic crossover time. This is consistent with the 
general discussion of a growing interface Q . Generally, 
the crossover time T co oc L z , where z is the dynamic ex- 
ponent. As in the linear case for A = discussed above 
the saturation width for a growing interface is a finite 
size effect and the transient growth does not saturate to 
stationary growth for an infinite system. Whereas the 
situation was easy to analyze in the linear case where we 
can identify the independent modes, and where the sys- 
tem size L is replaced by the inverse wave number 1/fc, 
i.e., the thermodynamic limit is probed in the infrared 
limit k — > 0, the situation is more subtle in the nonlinear 
soliton case since we dont have a normal mode structure 
but only approximate elementary excitations. 

On the other hand, for a finite size system, imposing 
for example periodic boundary conditions as indicated 
in Fig. 7 in order to ensure a growing interface in h as 
the slope solitons revolve, the velocity v becomes endowed 
with a scale and is quantized in units of L /T. Notice here 
the important difference between the diffusive case and 
the solitonic case. In the diffusive case the excitations 
are not propagating but are linear combinations of grow- 
ing and decaying modes as discussed in Sec. Ill, and the 
system size L only enters in the quantization of the wave 
number k oc 1/L, yielding the crossover time T~~ oc L 2 . 
In the solitonic case the localized modes are propagat- 



ing giving rise to genuine nonequilibrium growth. The 
system size L then enters together with the time T in 
setting a scale for the velocity v. 

A simple estimate, replacing the so liton amplitude 
u = u + — u_ in the general expression ( |5.30D by u cm — 
(u + +w_)/2 = —v/X from Eq. (5.26) and moreover v by 
L/T we obtain S so \ oc const. vh A j\ 2 T 2 which inserted in 
P oc exp[— S'soi/A] yields the soliton crossover time 



Tl 



sol 



L 3/2 



(5.31) 



and we infer the dynamic exponent z = 3/2. In the tran- 
sient short-time regime T <C T*° l the soliton configura- 
tions contribute to P; in the long-time regime T 3> T™ 1 
the soliton contribution vanishes and only the diffusive 
modes and their interaction s con tribute to P. We also 
notice that the expression ( [5.31 ) is consistent with the 
dimensionless argument \(A/v) i / 2 t/x 3 / 2 in the scaling 
function for the slope correlation function discussed in 
paper II and in QM^. In Fig. 8 we have depicted 
the crossover regimes. 



T diff , 2 
I co L 




FIG. 8. Here we depict the crossover time as a function 
of the system size L. In the early time regime for T <C T^" 1 
the distribution is dominated by soliton contributions. In the 
intermediate time regime for T c d ifi > T > Tco 1 the soliton 
contributions become suppressed leaving the diffusive mode 
contributions. Finally, for T 3> T^ S the diffusive modes also 
die out and we approach the stationary distribution. 

The short-time probability distribution P{h, T) for 
T <C T c s °' has been discussed within the directed poly- 
mer approach [ ^l|j66| . Based upon a replica scaling 
analysis [ |67| one finds for positive h (measured rela- 
tive to the growing mean height) the heuristic expression 
P(h,T) oc expl-Qi/T 1 / 3 )*], where n = 3/2; for h < 
based on numerical results, rj ~ 2.5. Recent exact results 
for the asymmetric exclusion model which falls in the 
same universality class as the Burgers equation, see ref. 
p8[ , where other references also can be found, also seems 
to have bearings on the height distribution. Using the 
Bethe ansatz method a skew distribution, characterized 
by the exponents 3/2 and 5/2, have been found for the 
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large deviation function of the time averaged current. 

Within the present soliton approach we can derive a 
qualitative expression for the early time height distribu- 
tion by noting that \u+ - u_| 3 ~ (ui) 3 / 2 (TA)~ 3 / 2 ~ 
/i 3 / 2 (TA)j 3 / 2 . Inserting this result in the general expres- 
sion ( 5.30| ) we obtain 



P st oc exp 



d d xdt p n 



dt 



(5.33) 



P so \(h,T) oc exp 



1 

AT 



1/2 



(5.32) 



in accordance with the directed polymer-replica based 
result and related to the exact results for the asymmet- 
ric exclusion model. The skewness of the distribution 
in h must then arise from the bias in the statistical 
weight exp[— S/A] assigned to the left and right hand 
solitons giving rise to a predominance of right hand soli- 
tons (S = 0), corresponding to relative forward growth. 
Unfortunately, our present understanding of the soliton 
approach and the inaccessibility of a more detailed multi- 
soliton solution do not allow a more detailed analysis. 



B. The general case 

In the general case for d > 1 the slope and noise fields 
u n and p n have lon gitu dinal vect or character and are 
governed by Eqs. ( [4.3j ) and (4.4) determining an or- 
bit (u n (x p ,t),p n (x pi t)) in p n u n phase space. At long 
times the orbit must pass close to the stationary sad- 
dle point (u n ,p n ) = (0,0) in order to induce Marko- 
vian behavior and then progress onto the stationary zero- 
energy manifold Eb = J d d x TLb = with energy density 
given by Eq. (4.2), yielding the weak no ise distribution 
P oc exp[— Sb/AP with Sb given by Eq. (O). 

In d = 1, as discussed in Sec. A, the orbit on the sta- 
tionary manifold is attracted to the submanifold p = 2vu, 
yielding the symmetric stationary distribution (|l.9| ) with 
time-dependent skew corrections. For d > 1 this behav- 
ior is only encountered in the linear Edwards- Wilkinson 
case for A = 0; note also the general discussion in Sec. 
II. The attraction to the submanifold p n — 2vu n is as- 
sociated with the underlying fluctuation-dissipation the- 
orem and gives rise to the stationary Gaussian distribu- 
tion P st oc exp[— (v/A) J d d x u n (x) 2 }; the corresponding 
free energy is F = (1/2) J d d x u n (x) 2 . This distribu- 
tion yields the roughness exponent £ = (2 — d)/2. The 
dynamic exponent z = 2, corresponding to the diffusive 
mode contribution; note that the Galilean invariance is 
not operative in the linear case and that we consequently 
dont have the scaling law constraint ( + z = 2. 

In the nonlinear case A ^ for d > 2 the long-time 
orbit emerging from the vicinity of the stationary point 
(u„,p n ) = (0, 0) diverges for larger u n from the subman- 
ifold p n = 2vu n which constitutes a sort of tangent plane 
to the zero-energy surface at the stationary point. In 
the limit A — > the energy surface then collapses to the 
tangent plane. More detailed, in the stationary limit the 
distribution is given by 



Noting that p n is slaved to u n on the orbit we obtain in 
general (to leading order in V„) 



Pn = u n F d (u 2 ) + u p V n u p Gd(u 2 ) 



(5.34) 



where the scalar functions Fd and Gd depend on the in- 
variant u 2 = u n u n and parametrically on the dimension 
d. In one dimension we have F\ = 2v and G\ — 0. For 
d > 1 to leading order in V n the distribution P st is even 
in u n - Note, however, that for Gd ^ there is a skew 
correction to P st - Assuming that this analysis is valid, 
the determination of Fd and Gd is the n give n by the so- 
lution of the equations of motion (4.3) and (4.4) on the 
stationary zero-energy manifold. 

In the nonlinear case for d > 2 the rough phase gov- 
erned by the strong coupling fixed point only appears 
for a renormalized coupling strength A = (A A) 2 /V 3 ex- 
ceeding a finite threshold value A c J69|-[7l[ and recent 
work |55|,[56| moreover indicates that unlike the case in 
d = 1 even the stationary probability distribution ex- 
hibits skewness, i.e., P s t{u) ^ P st (—u). Such a behavior 
does not seem compatible with the analysis above and 
indicates that the present weak noise approach only ap- 
plies to the weak coupling phase for d > 2. Presumably, 
the strong coupling phase for A > A c is only accessed be- 
yond a critical noise strength A c , A c = (A C A) 2 /V 3 . These 
issues and the fact that the strong coupling fixed point 
(fortuitously) can be analysed in d = 1 in a nonpertur- 
bative weak noise approximation remains not very well 
understood and call for further investigations. 



C. Summary and conclusion 

In the present paper we have advanced a general weak 
noise canonical phase space approach to stochastic sys- 
tems governed by a Langevin equation driven by addi- 
tive white noise. Reformulating the associated Fokker- 
Planck equation in the nonperturbative weak noise limit 
in terms of a Hamilton- Jacobi equation we have discussed 
the time-dependent and stationary probability distribu- 
tions from a canonical phase space point of view. The 
issue of solving the stochastic Langevin equation or the 
associated Fokkcr-Planck equation is thus replaced by 
solving coupled Hamilton equations of motion determin- 
ing the orbits in phase space. The stochastic nature of 
the underlying problem is reflected in a peculiar topology 
of the energy surfaces different from the one encountered 
in ordinary dynamical problems. The Markovian behav- 
ior thus corresponds to the existence of a stationary hy- 
perbolic saddle point which controls the behavior of the 
orbits in the long time limit. 

We have in particular applied the canonical phase 
space approach to the noisy Burgers equation describing 
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a growing interface in one dimension. We have recov- 
ered the well-known stationary distribution and derived 
expressions for the time-dependent distribution, at long 
times governed by linear diffusive modes and their in- 
teraction and at shorter time by nonlinear soliton exci- 
tations. In higher dimensions where the noisy Burgers 
equation predicts a kinetic phase transition to a strong 
coupling phase the canonical phase space approach only 
seems to access the weak coupling phase. 
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